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ABSTRACT 

A  mathematical  model  of  a  branching  stochastic  process  utilizing 
generating  functions  is  presented„  The  probability  distribution 
of  the  number  of  members  of  the  process  at  discrete  time  periods, 
z^j,  and  the  probability  of  extinction  is  discussed^  When  there  is 
a  non-zero  probability  of  surviving  indefinitely,  the  norraed  random 
variables,  Zjj/E(zn),  converge  with  probability  one;  the  cumulative 
distribution  of  this  random  variable  is  absolutely  continuous „  The 
time  until  extinction  and  the  total  number  of  members  of  the  process 
is  examined  when  the  probability  of  extinction  is  one^  The  distri- 
bution of  the  z^  given  that  z^^  is  not  zero  is  discussed  for  this  case„ 
The  maximum  likelihood  estimates  for  the  probabilities  involved  in 
the  process  are  determinedo  An  example  is  given  of  a  branching  pro- 
cess  in  which  the  probabilities  are  dependent  on  time  and  a  solution 
is  found  using  Laplace  transform  methods. 
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NOMENCUTURE 

gymbol       Meaning 

P(A)         the  probability  of  the  event  A 

P(A  B)        the  conditional  probability  of  A  given  B 

P(A,B)        the  joint  probability  of  A  and  B 

Zj.  the  number  of  members  of  the  r^^  generation 

X^^"^'         the  number  of  offspring  of  the  i   member  of  the 

nth  generation 

the  expected  value  of  Zj^ 

the  variance  of  z, 
f(s)         generating  function  of  z-, 

Pj.  probability  of  r  members  in  the  first  generation 

fn(s)        generating  function  of  Zj^ 

Pj^          probability  of  r  members  in  the  n''^^  generation 
a           probability  of  extinction 
w„  a  normed  random  variable  ( £i 


n 
w  the  limit  of  w^^  as  n  approaches  infinity 

Gjj(u)  the  cumulative  probability  function  of  w^^ 

j^(s)  the  characteristic  function  of  w^ 

G(u)  as3nDQptotic  distrij^ution  branching  from  f(s) 

k(s)   .;''  probability  generating  function 

H(u)  the  asymptctic  distribution  branching  from  k(s) 

(s)  the  characteristic  function  of  H(u) 

a.^  probability  that  z^  j^  0 

p^  probability  that  z^  equals  r  given  that  it  is  not 

equal  to  zero 

V 


jTjj       number  of  members  of  the  n^^^  generation  given  that  it  ij 

not  zero 
Fj-^(s)     probability  generating  function  of  j^^ 
N        number  of  generations  to  extinction 

(s)      moment  generating  function  of  N 
S        total  number  of  individuals  in  the  process 
g(s)      probability  generating  function  of  S 
Zj^j.      number  of  individuals  in  the  m^^  generation  who  have 

exactly  r  descendents  in  the  (m  1)   generation 
Sjj       total  number  of  individuals  in  the  first  n  generations 


vi 


INTRODiJCTION 

The  present  paper  is  an  essay  on  the  theory  and  application  of 
discrete  branching  stochastic  processeSo  Most  of  the  tsyii^totic 
theory  was  extracted  from  the  work  of  Harris  [lO]  who  derived  many 
of  the  relationships  for  this  class  of  stochastic  processeSo   Let 
us  introduce  these  processes  by  means  of  a  historical  exairpleo 

In  l67li,  one  of  the  first  investigations  into  discrete  branching 
processes  was  made  by  the  Rev.  H.  W,  Watsono  Dr„  Fo  Galton  had 
posed  the  problem  of  extinction  of  family  names  to  the  mathematical 
community  in  England  since  many  names,  prominent  in  history,  were 
no  longer  to  be  found.  The  problem  was  posed  as  follows:  Assume 
a  large  nation  with  N  adult  males,  each  with  a  different  surname,, 
Wliat  proportion  of  names  would  be  extir^ct  after  r  generations?  In 
the  r'th  generation,  how  many  people  would  have  the  same  name?  Rev„ 
Watson  [8]  derived  some  relationships  that  enabled  him  to  coirpute 
the  probabality  of  extinction  and  made  some  sample  calculations  to 
demonstrate  the  difficulty  of  computing  numerical  answers^   Notable 
among  these  relationships  was  the  functional  relation  of  the  proba- 
bility generating  functions  involved  (section  III,  theorem  1) o 

Branching  processes  received  some  attention  thereafter  and  the 
theory  was  applied  to  problems  of  gene  mutation  [7],  population 
growth  [12]  ,  cosmic  radiatioii  [l]  ,  ionization  fluctuations  [l8l , 
neutron  production  in  atomic  piles  [llj  ,  and  cascade  theory  [20]  o 

In  this  paper  we  will  review  the  theory  of  simple  discrete 
branching  processes  and  some  of  its  applicationSo  Let  us  consider 
the  original  problem  of  survival  of  family  names  and  state  the 
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basic  problem  of  branching  processes  in  that  frameworko   Suppose 
an  immigrant  by  the  name  of  Door  comes  to  the  United  States  and 
there  is  no  one  else  here  having  that  surname,,  What  is  the  prob- 
ability that  the  name  survives  r  generations  and  how  many  men  vdll 
there  be  having  that  last  name? 

There  is  a  probability  that  Mr^  Door  will  either  not  marry  or 
not  have  any  male  children;  there  are  other  probabilities  associated 
with  the  number  of  sons  he  may  have„  We  will  denote  these  probabili- 
ties by  Pj,,  i.eo,  the  probability  of  no  sons  will  be  Po,  of  r  sons, 
Pj,,  ViJe  need  a  symbol  for  the  number  of  individuals  with  the  name 
of  Door  iii  the  n"^^  generation,  so  let  z^^  be  the  number  of  individuals 
with  the  name  of  Door  in  the  n"th  generatioHo  Then  we  can  write 
P(Z]_  2  r)  =  pj.  indicating  that  the  probability  of  r  individuals  with 
the  name  of  Door  in  the  first  generation  is  pj..,   If  Mro  Door  coes 
have  sons,  each  of  them  may  have  children  and  tte  process  is  continuedc 
This  may  be  represented  as  a  tree  type  graph  with  an  initial  node 
representing  his  sons,  and  so  ono  A  method  of  analysis  of  the 
problem  based  on  graph  theory  was  proposed  by  Otter  [19]  o   In  this 
paper  the  method  of  generating  functions  is  usedo 

We  will  assume  that  the  number  of  sons  born  to  one  of  Mro 
Door's  sons  follows  the  same  probability  law  as  that  of  Mr,,  Door., 
In  other  words,  the  probability  that  Door,  jTo,  has  n  children  is 
Pj^o  Then,  if  Xj^^^^  is  the  number  of  offspring  of  the  i^^   individual 
in  the  n^-h  generation,  we  have  P(zn+i  =  r|  Zn  s  m)  s  P(^  x]^"'  =  n)  „ 

This  is  the  outline  of  the  problem,  the  notation,  and  the  assuiip- 
tions  made  in  this  paper.,  Before  taking  up  the  mathematical  theory, 
some  other  problems  will  be  posed  in  this  frameworko 
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In  atomic  piles,  the  number  of  free  neutrons  is  an  imjoortant 
factor o   If  we  start  with  one  neutron  and  a  fission  takes  place,  the 
probability  that  r  neutrons  are  released  is  a  nuclear  constant  which 
must  be  empirically  determined  for  the  amount  and  arrangement  of  active 
material.  The  probability  of  no  fission  (no  descendents)  is  the 

average  probability  of  leakage  or  absorptiono  An  atomic  pile  in 

« 

a  nuclear  reactor  is  subcritical  when  the  reactor  is  not  operating 
and  there  is  no  possibility  of  an  explosion,  ioBo,  the  probability 
that  the  number  of  fissions  increases  sufficiently  to  cause  an 
explosion  is  zero;  the  mean  of  the  probability  law  for  a  subcritical 
mass  is  less  than  one.  As  the  amount  of  fissionable  material  is 
increased,  the  probability  of  an  explosion  finally  becomes  greater 
than  zero.  liJhen  this  occurs,  the  mass  is  termed  supercritical  and 
the  mean  of  the  probability  law  for  this  situation  is  found  to  be 
greater  than  one,  A  nuclear  reactor  is  termed  critical  Vhen  it  is  , 
operating;  the  mean  of  its  probability  distribution  is  exactly  oneo 

The  organization  of  the  paper  is  as  follows:  The  general  math- 
ematical character  of  a  discrete  branching  process  is  presented  in 
section  lo   It  will  be  seen  that  the  mean  number  of  progeriy  per 
individual  plays  an  important  ro3,e  in  the  asymptotic  theoryo  This 
theory  is  presented  in  section  II  for  the  mean  of  the  probability 
distribution  greater  than  one  and  in  section  III  if  it  is  'not 
greater  than  one.  The  maximum  likelihood  estimation  of  the  distri- 
bution  of  progeny  per  individual  is  given  in  section  IV.  An  exairple 
of  an  ionization  cascade  is  presented  in  section  V.  This  is  a  more 
general  branching  process  than  those  treated  in  the  model  used  in  the 
earlier  sections  and  its  treatment  suggests  how  the  methods  of  Markov 

processes  may  be  applied. 

-3- 


GENERAL 
The  structure  of  branching  processes  was  presented  in  the 
introduction  where  several  assumptions  were  made  aiid  some  notation 
was  introduced.   It  is  appropriate  uo\-i   to  summarize  the  assumptions 
made  in  this  paper  and  to  extend  the  notation.. 

Consider  two  related  classes  of  random  variables: 

Zjy       the  number  of  "individuals"  in  the  n^^  generation 
X^^^   the  number  of  "sons"  of  the  i^^  individual  in  the 
nth  generation, 
They  are  related  by  the  equations : 

Zm  -  A  Ai  for  n=  1,2,3,.  o  o 

The  quantity  Zq  is  the  number  of  the  initial  group  who  start  the 
process;  Zq  is  assumed  to  be  one,  since  this  does  uol   affect  the 
generality  of  the  theory  as  is  pointed  out  at  the  end  of  this  sectiono 
The  probability  distribution  of  the  first  generation  is  given  by  pj,, 
ioe.,  the  probability  that  there  are  r  members  of  the  first  generation 
is  pj,. 

The  following  assumptions  are  made  throughout  this  paper: 
(1)  Given  Zj-^  is  equal  to  m,  the  random  variable,  Znfi  ,  is  the  sum 
of  m  independent  random  variables,  each  one  of  which  has  the  same 


distribution  as  z^  :   Pi^L     -  '"/  ~  Pr 


(2)  The  variance  of  z-|_  is  finite. 

(3)  Pj.^  1,  for  all  r  j  otlierwise  the  process  degenerates  to  a  simple 
arithmetic  calculation.  ; 

(h)   Pq+  P]_  ^  1,  otherwise  we  have  the  binomial  case  whose  treatment 
in  the  context  of  branching  processes  is  trivialo 


The  following  notation  is  used  throughout: 
y^^E-i^i)-;^  rp^  ^      the  mean  of  zi 

-j_  CO 

(T^^Vori?,)  =^(r^p^j-^a    the  variance  of  z^ 
ao 
~  ^0  ^^^^  probability  generating  function 

of  Z]^  (s  is  a  complex  variable)* 

oo 
"^  Pnr"^  "^^^  probability  generating  function 

of  z^o 

From  the  assiimptions  listed  above,  it  is  clear  that  : 
f  (s)  =  s,  p,   =  p  ,  f  (s)  and  f  (s)  are  continuous  on  the  set  of 
points  consisting  of  | s | < 1  and  s=lo         ' 

Several  authors  (see  [2]  ,  [9]  ,  p-lj  )  have  examined  the  generating 
functions  and  have  found  some  basic  relationships,   Some  of  these  are 
presented  below. 
Theorem  1^  fi-,(s)  is  the  n'^'^-  functional  iterate  of  f  (s)  : 

fn(s)  =  f(Fn.,(s))  =  Fn.|{F(siJ 

Proof:  From  the  irdtial  assumptions,  the  X^  '   are  independent  for  all 

It  follows  that  I  (^n+l"'^i^n-J)      is  the  coefficient  of  s^ 

in  the  power  series  expansion  of  (f(s))"^  since  the  generating  function 
of  the  sum  of  j  mutually  independent  random  variables  is  the  product 
of  their  generating  functions.   See  Felier  [6,  p.,  250ff ,]  , 
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Hence,  p  (s)  =  ^  p  .  (flslf  =f(f(Sl) 

oo 


fnls)-2  Pn-i^f^5^^^  ""  ^n-it^fS)]  -F,.,lfL.  .  (f(f(S)))  ...)] 

The  coefficient  of  s^  in  fj-^(s)  is  the  probability  that  there 
are  r  individuals  in  the  n''"  generation,,  Thus  the  questions  posed 
by  Dro  Galton  ma;/  be  answered  in  the  framework  of  the  present  models 
The  constant  term,  V^o^   ^^   ^^®  probability  that  the  family  name  is 
extinct  by  the  n'^^  generation. 

The  mean  and  the  variance  of  z^^   can  be  determined  using  theorem  lo 
Theorem  2 o   p-  /..  i    «      w   /,  \«  S^  iuJ'iU^''\\ 


-  r\6^  /^- 


E^=    fr;(i)  =  |rp^^  =E12.J 


T'O 


Therefore,     Var  (2n)^  'r^  (')  ^  't^U)  -  [f^^'iOj    . 
By  theorem  1         Fk,^,  (  S)  -  F^  Cf  (s)]    -ffF^ISJJ 

Taking  derivatives,      f^^,(s)  =  f^  CHsljf  (S)    -  f '  [f^  (s)]  f^  (S) 


=  f"CfJs)3Cf„'(s)]^    +  f'[F„(s)]f;(SJ 

Hence         Ffm  I '^  " /^  ^  ( ''"  1   "  M 


Solving  for  f{.^  (1)   gives 


•^^.i*^  /  ivh 


"Thus  Var ( z.J  r:     ^  /^"^  (^/U"-') for    ^  ^ ^ 


For  the  case  yU*!  ,    L'Hopital's  Riile  can  be  used; 


3./U  - 

and     Var   ( z  J  =    n  5*^  :f or     yU  - 1 


Higher  moments,  if  they  exist,  may  be  found  by  a  similar  pro- 
cess utilizing  the  higher  derivatives  of  ^nfl^^^" 

A  subject  of  great  interest  in  any  branching  process  is  the 
probability  of  extinction  of  the  specieSp v  The  Fundamental  Theorem 
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of  branching  processes  tells  us  that  this  probability  is  a  fixed 
point  of  the  generating  function.  The  proof  given  follows  that  of 
Bharucha-Reid  l?]"     First,  we  ixote  that  if  "a"  is  the  probability  of 
extinction,  then  aslim  p^-,s  lim  f„(0)o 

Theorem  3.   The  probability  of  extinction,  a,  is  the  smallest  pos- 
itive number,  p,  such  that  : 


Further,  a  s  1  if  and  only  if  yU  ^  1, 

Proof ;  We  see  ^that  a  =  0  if  and  only  if  p^  ~  0  and  f (1)  =  1,  so 

that  0-Q  -1  as  it  should  be^ 

To  show  that  a  =  f(a),  we  note  that  since  p^.^^  ~   fj^(O), 

f\^90     no      lo-ioo     ^^  n-»oo         n-i  n-^^o         n-XyU 


so  a  = 


Now  f (s)  is  a  power  series  in  s  with  positive  coefficients  and  hence 
strictly  increasing  for  real  arguments,.   If  0<a<b,  then  f(a)^f(b)o 
Also,  since  Pno^^  ^^^  ^1^  ^^5  then  a2  0. 
Let  p>0  be  such  that  p=  f(p)<>   It  follows  that 

Assumirig  P^o- P>  then  fn(0)^p  and 

So,  by  induction,  p^o'^P  ^^r  all  n„ 

Then  lim  Pno-P  ^^^^   a^p. 

It  remains  to  show  that  a  =  1  i"f  and  only  if  yH   ^  j  ^ 


Assume  first  that  a  =  1  and  show  that  u   ^  \ 
Now,  it  may  be  noted  thai : 

f  (s)2  0,  f"(s)2:0,  so  f(s)  is  convex  in  (0,1) 

f(0)  s  Pq>0  ,  f(l)  =  1 

f(s)y»s  for  s  in  (0,1)  since  1  isthe  minimum  positive  number 
satisfying  s  =  f(s). 

Therefore,  1  -  f(s)^l-s  and. 

■.°  ■  ■    is  bounded  by  1  and  is 

monotonically  increasirig 
with  s  for  s  in  (0,1). 
Therefore,  f'(l)  exists,  f'(l)^l  and  ^^i. 

Now  assume  ^<lo   Then,  since  f(s)  is  convex,  f'(s)  is  either 
constant  in  (0,1)  or  strictly  increasing  with  So  By  assumption, 
p.^  ^  1  so  p  >0  and  f  (s)  ^  1  in  either  case  for  s  in  (0,1)  „ 
By  the  mean  value  theorem,  ^  K\\  '  ^[^)-  f  (X\(h^)       wkei^e  s^^^i 


Then,  /-F(s)  =  r-'()r)(|-S) 
F(S)>S     for 


and 


O^S^i 


Therefore,  a  =  1. 

Using  theorem  3,  it  is  a  simple  matter  to  determine  the  limits 
of  the  Pj^j,  for  r  5^  0„ 

Theorem  Uo  A  discrete  branching  process  either  dies  out  or  becomes 
extremely  large,  i.eo,  lim  p.^^  ~  0   for  r  ^  Oo 


Proof:  Let  t  be  a  real  number  such  that  Oi' t<  1  and  recall  the 
iteration  property  of  fn(t) » 

Case   1.^41      Then  f  (t) ;?  t   (see  proof  of  theorem  3),   and 
f(t).^f2(t)^  „    „    „   ^fn(^)^  o    o    .    <1. 
Therefore,   lim  f^ (t)   s  1  ^  a. 

Case  II,^>i  and  t^a  (Pq>0)o     Then  f(t)>t  since  f(s)   is 
convex  and  f(t)^  f2(tK   »    •    <>    -^^n^^^^  °    °    ' 
Therefore,   lim  f  (t)   =  a. 

Case  IIIoyu.>|  and  t>  a. 
Then  f(t)4t  due  to  the  convexity  of  f(s)„ 
Hence  f(t)>f2(t)>.    ,    „    >fn(t)>   o    o    o  >  a„      p 
Therefore,   lim  f„(t)   =  a. 

Now  fn(t)  =  ^  p^^t      -  p^^  f  2  f^r^         ^^  definition.   Taking 

the  limit  of  the  last  expression  gives j 


lim  f^(t)   r  lim  p      +  lim  >>  Pv^v.'t'''     =  a  f r  .m  the  above  three  caseso 

OO 

Therefore,  lim  ^>  p„„t^  =  Oo   Since  the  series  converges  uniformly 
in  an  interval  one  may  pass  to  the  limit  before  summingo 


oo 


Therefore    ">  lim  Pr.-rt^  ^   0  and  lim  p^,_  =  0  for  r  ^  0  by  the 

uniqueness  of  power  series. 

Theorems  3  and  U  give  us  some  interesting  information^,   If  we 
are  working  with  an  atomic  pile,  the  probability  is  one  that  it  will 
die  out  unless  it  is  suoercritical;  then  tlie  probability  that  it 
will  blow  up  is  1  "  a„ 
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Lotka  [17]  found  that  the  probability  law  of  the  number  of  sons 
a  man  might  have  was  closely  approximated  by  a  type  of  geometric 
distribution,  namely,  p^  =  OoU825,  Pj^  -  (Oo2126)  (Oo5893)^"-'-   klTl 
based  on  a  1920  census „  '  It  follows  that: 


f(s)  ^  0„1|825+  Oi2126s     ^^  f«(s) 
l-0<,r"93s    ' 


=  lo2622 
s=l 


Solving  f (s)  =  s,  it  is  found  that  a  s  0o6226o  Based  on  these  values 
the  probability  of  extinction  of  a  family  name  is  (0o6226)^^  where 
n  is  the  number  of  male  members  of  the  family  at  the  present  time.. 
This  formula  follows  from  the  fact  tha\  if  there  are  r  members  of 
the  zero"^*^  generation  then  the  generating  function  for  the  first 
generation  is  (f(s))r.  Then  fj^(s)  *  fn_i((f (s))^  ?  if  f(a)  =a,  then 
(f(a))^  s   a^  and  a^  is  the  probability  of  extinction. 
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MEAN  GREATER  THAN  ONE 

In  this  section  we  will  consider  the  case  where  >^  is  strictly 
greater  than  one  and  the  family  has  a  non-zero  probability  of 
surviving  indefinitelye   Since  we  have  determined  that  z-^^   may 
become  extremely  large,  we  will  work  with  a  normed  random  variable^ 

Zn 
Let  Wj^  ~   — Tj-  o  The  mean  and  variance  of  Wj^  are  easily  deter- 

minedo 


£(-^)-E(#„)=>r-^~w^^y"i 


Treatment  of  the  convergence  of  w^  may  be  found  in  Harris  [lO] 
and  Bharucha-Reid  [2^  o     We  will  prove  first,   the  convergence  in 
probability  and  then  the  convergence  with  probability  one^ 
Theorem  5o      IfyU>l,   the  random  variables,   w^,    converge  in  probability 
to  a  random  variable  w. 
Proof o      Let  n,   m,   be  integers  and  n^-mo 


E  ("ZnlZnj-r)  =  X  J  [coefficient  of  sJ   in(fn.in(s))^] 


J=l 


.12- 


5  =  1 


r?o    "^i^ 


/u"-'"E(2;S) 


£(w/' 


ho 


E[(u>n-u>^)^]  ^    Etu;')- E(u;.^i 


w^  h-»o© 


■'^     Eju;^-l/>mfj   ^  0 


Hence  the  random  variables  converge  in  mean  square  and  by  a  theorem 
of  Kolmogorov  [^13,  p.3U,    l]  ,   in  probabilityo 

Convergence  with  probability  one  may  be  proven  with  the  aid  of 
Doob's  [U]  Martingale  convergence  theorem^ 


.13- 


A  discrete  parameter  martingale  is  a  stochastic  process  such  thats 
^)  2  dyj)  <   «3    for  all  n 

^'^  ^  ^^n4l'  ^n'  yn-1,  »  o  o  ^  Jl)  ^  /n 
The  theorem  may  be  stated  as  followSo  Let  {y^^   nZQ)   be  a  discrete 

parameter  martirxgale^  Then  E(  jyj  )^  E(|yJ  )  i  E(  72  )4  o    o    o 

If  lim  E(|y„|)  =  k  ^''O   then  lim  y^  s  y*  exixts  with  probability  one 

and  E(|y*|)  ^  ko 

Considering  our  random  variable  w^  ^  E(wn)  ^  1  for  all  n, 


-Wr 


Hence  the  (w^)  form  a  discref  parameter  martingale  and  we  have 

Theorem  Uao  If  iA>lf   the  random  variables,  w^^  ,  converge  with  probability 

one  to  a  random  variable  Wo 

Since  Wj^  converges  to  w,  we  may  discuss  the  probability  distri- 
bution  of  Wo 

Let  Gj,(u)  -  P(wniu),(pj,(s)  ^  E(  e"n^)  ^  J^e^^  dG^(u), 
G(u)  s  P(w£u),  lp(s)  s  ]  e^^  dG(u).  The  distribution  G(u)  can  be 
called  the  asymptotic  distribution  branching  from  f(s)o  There  is  a 
very  interesting  functional  relation  between  the  moment  generating 
function  fP(s)  of  w  and  the  probability  generating  function  f(s)  as 
given  by  the  following  theoremo 


aU- 


Theorem  $.  L()  (/XS)   =  F(Lp(S)) 


Itmn    Gn^^)  "  G  t^)   .therefore, 
h-»oo 


Upon  taking  the  limit  of  (f^^,  (yUS)=  '^TlP^  (S)]   the  desired  result 
is  obtained. 

It  also  follows  that  the  k^^  moment  of  w  exists  and  is  finite 
if  and  only  if  the  k"^^  moment  of  z-j_  exists  and  is  finiteo 

Since  G„(0)  ~   P(wni  0)  s  p     then  lim  0^(0)  "  a  ard  P(w  s  O)sao 
The  following  theorem  will  be  used  later „ 

Theorem  6o  Let  Gn(u),  G2(u)  be  distributions  having  equal  first 
moments  and  finite  second  moments  such  that  their  characteristic 
functions  ^^( it)  and  (^2^^^^  satisfy  Li?j.  (tt^l^  P[^j.(^^ij  r  =  1,2 
Then  G-j_(u)  =  G2(u) 


Proof:       if)  {it)   ^\    ^     ^     ±L    (ttj^   0[i''] 

for    }tj    small  by  a  special  form  of  NcLaurin's  theorem  for  small 
values  of    !tL      See  [Cramer     3,  Po27J    o     Then, 


i  a?,(^ti-  a);i(a)|  =  11+  J;  ji  (ar  f  o{t^)  =1 


a       .UJ 


^  *o 


-5"^(Urt0(£')     where    /^/^'(lt)  =  lV 

Therefore,     |  L^,  itt)  -  ifg  (  U  )  |  ~  t'^/S  (tj        \jj^zrt     \\m     /Sii) 
Then      I  0^,  (li/>L)  »  LP^  (/^Uj)  =  |  F  Ctfj  l^i)J-  FLCgilUiJi    by  theorem  ^ 

^^  lll'jI.L^i'l'Sl^i*^!       sincefrts)!^^  when  Isl^J 
Hence  j  t'^/A^ S  (/Ut )  j  £  yU  I  ^''^U)  i  .  •   .. 

and       yu|(d(/Xt)j  1^   !|S(t)|  implying  that      >U  i/<?(fl)N    I  (3lo)\ 

But JLA>1,  SO   ^(t)  must  be  identically  zeroo 
Therefore     |  l^j  (i.t) -Cpa.!^^)!  =  0 

and       (P,  (tt)i  LPj^llt) 

Finally  G]^(u)  s  G2(u)  by  the  uniqueness  of  moment  generating  functionso 

In  order  to  determine  some  properties  of  G(u),  it  \id.ll  be  con- 
venient to  work  with  a  closely  related  function  H(u)  which  is  defined 
as  the  asymptotic  distribution  branching  from  k(s)  where 


To  demonstrate  that  k(s)  is  a  generating  function  use  the 
fact  tliatyu>l  implies  f(a)  =  a^  1.  Then 

I  —  CK 

Clearly  k(s)  is  monotonically  increasing  since^  if  t ^  s,  then 

Therefore  k(s)  is  a  probability  generating  function^ 

It  is  interesting  to  note  the  following  properties^ 

H'ii)   =  r'(!)  ^^ 

K'  (S!=  r'[s(i.a)fa](!-a] 

The  coefficient  of  s^  is  (l^a)*^~-^(p  ^  f  ^  'Ka  -f-  .  .  o  f  |  |^  j  P   C« 


Q     f  .  .  . 


Let  lj;(S)  s  f  ^*^'JH'»uJ   °  ^^  ^  ^^^^®  °^   argument  similar  to 
the  proof  of  theorem  5^  it  can  be  shown  that  {^  (jU-S)   =  K  L  ipi^)j. 

=17" 


To  demonstrate  that  k(s)  is  a  generating  function  use  the 
fact  thatyU>l  implies  f(a)  :  a4  !<,  Then 


I  -a 


f</|j  =  iilLi^L^^  i 


-a. 


Clearly  k(s)  is  monotordcally  increasing  since j,  if  t-<s,  then 

/^^  1  —  ol 

Therefore  k(s)   is  a  probability  generating  functionc 

It  is  interesting  to  note  the  following  properties c 

The  coefficient  of  s^  is   (l-a)k-l(p^  f  p^^^^Ka  f-  .  ,  -  f  ("'^1  fir.h''''^ 

Let     l|^(S)    -    f   J^^^(^'^V)      °   ^^  ^  ^^^^®  °^  argument  similar  to 
the  proof  of  theorem  5^   it  can  be  shown  that     W  iU-Sf  =  K  L  y^iS)]. 

=,17„ 


We  now  proceed  to  determine  a  representation  for  the  random 
variable  Wo  First  we  determine  the  relationship  between  l^(s)  and 
(pis). 
Theorem  7.      ^(S)  -  ^^-77^-"- 

Proof:  Let  ijJjiS)^      (PC^-CcJS]  -  g. 

I  -a. 

Then  -  (|/;  (S)  =  MlUtM}'1i^'0,)_    ^    (p'ro-ajsl 

and   i/^/(^i~  ^^io)-^  ^     Therefore  i^s   and  ^'  have  equal 

first  moments.   Since  ^^"(sj     s  (^'^[d-als]  (s-a)  j    <  o^ 

Ij^i  and  Lj/'  have  finite  second  morentSo 

Therefore  y  and  \^   satisfy  the  same  functional  relationship  with  k(s)o 
Therefore  by  theorem  6,  i^  (S)  s  tp^j  (s). 

Now  we  can  determine  the  relationship  between  H(u)  and  G(u)c 
Theorem  6.    Hiv)  ^    6  (l^g  ]'  0.  ' 

Proofs  From  theorem  7,   i|^  (S)  r  '^pPHQ-JSj-a 

/-a 

By  definition  of  generating  functions   ||/?(S)  -  J  t^^ d  Hlu) 


and  ...i^j  j'°°e"'''"'''*dS{u)  -a  j  =  ^L('t:a'sjj,a 

Hence  I 


Jo 
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Let  EM  be  defined  as  follows:         £  f  ul  s  0        f  Of      U  ^  0 
Then       ^=J        d  2(u)    ^    j^M^^dc(u), 


Therefore  a  is  the  Laplace  transform  of    fivl. 


Hence    f  f^"  ( l-a)d  W(W   =    f  ^  "'^  ^  ^  i^*)  "  i  ^ '"'^"'^^ 

Therefore    ||-a)  H(.u)  =  S  f  pa)  "^  ^'^-^ 

and       H(l^)  =  -yia'  [S(l'^|~^]        for    U^-O 
Clearly  H(u)   s  0  and  H(«» )   =  1  as  it  shouldo 
To  put  the  relationship  of  H(u)   and  G(u)   in  terms  of  random 
variables,  we  prove  the  following  theoremo 

Theorem  9o     The  random  variable  w  is  distributed  as  the  product  of 
two  indepei  dent  random  variables,  w^  and  w',   where  P(Wq  s  q)   ~  a, 
P(Wq  s  _± — ,)   s  l_a  and  w»  has    &he  distribution  H(u)o 

Proof:    £(£^'^'sj  s  a  +(i-a)^"e.-i  J  HW 

=  a  ^■ll-ajC"e"<JH('/(l-(l)) 

0 

f 

-  tp(s) 


=19^ 


So,  P(  ujii-a)  ^u\uj>o]-=  H(u). 

For  later  use  a  theorem  due  to  Harris  [lo]  is  given  with  some 
remarks.,  Let  J=  lo^  1_)  ^  lo|a(~^j)<.   If  q-^  ^  0,  let  /=  oo  . 

Theorem  10,  If  ar<^,  Re{s)  £  0   ,  s  ^  0,  then  U)  {s}   =  ^^^   -f  Mo(s) 

where  M(s)  is  continuous  for  s  j^  0,  KiUs)   =  M(s);  Mq(s)  s  0  ("-^y j  iSl-*"^ 
M(-s)  =  M(s)o 

Remarks:  (a)  Under  the  conditions  of  the  theorem  M(s)  is  real  and 
poaitive  when  s  is  reg.1  and  negative „ 

(b)  If  ^{l^j'ioo,   the  r^'^  derivative  of  Ij^(s)  satisfies 

(c)  If  y=  oo  ,  IV  (s)  and  as  maiiy  derivatives  as  exist  approach 
zero  exponentially  as  |s|-»c70o 

We  will  now  determine  some  of  the  properties  of  G(u) ;  however, 
it  will  be  convenient  to  prove  theorems  for  H(u)  and  then  interpret 
the  results  in  terms  of  G(u),  Let  h(u)  =  H  (u) „ 
Theorem  11,  H(u)  is  absolutely  continuous. 

Proof ;  Cramer  [3]  breaks  up  an  arbitrary  cumulative  distribution 
function  F(x)  into  three  pieces:  .  /. . ' 

Fix)  =  Q|  F|  (x)  f  0^  F;j.()f)  -f  0^^  Fj  Ix) 

where   O ,  +  C«^  f  CI3  =  J       ^i^O  j     1=1^2^3 
F^ (x)  is  absolutely  continuous, 
FpCx)  is  a  step-function  and  is  equal  to  the  sum  of  the  saltuses 

of  F(x)  at  all  discontinuities  which  are  less  than  or 

equal  to  x, 

F.(x),  the  "singular"  component,  is  a  continuous  function  having, 

almost  everyvhere,  a  derivative  equal  to  zero, 

-20- 


LiXy^ 


O^  IP2  ! -J  ~  y-''    p   6  is  the  sum  of  an  absosutely  coiivergent 

rJVo  "^ 

trigonometric  series  and  is  thus  an  almost  periodic  function  which 
comes  as  close  to  b.^^   as  we  please  for  arbitrarily  large  values  of  t< 
By  theorem  10,  lim  tP  (t)  s  Oo  Therefore,  for  H(u),  a^  s  0  and  H(u) 


is  continuous o 


dv  ^  e 


^>^i,^j"  ' 


Let 


h^^(u)  = 


I.   o"^^^ 

V  =  ~i- e 


1;=  {j^iit) 


-im 


'W 


I 


Liu 


^TTLU 


Ji 


(A) 


^.M=-■5F^^[^■Me■'""-^'-^-'^'i+f 


e 


^{phj) 


"'-KM  "  * 


Examine  each  term  of   (A)   as  m  -*  00 


First      '\j^iL/n)e 


^Imiu 


MUim) 


'rni 


p-    -t   /lo(s)    e"'*^*^       by  theorem 


10  0 


.-y 


Imf    IMim]  i   Of- 


fe.)  1 


-Lini^ 


As  we  pass  to  the  limit,  this  term  is  a  function  of  continuous 
functions^  hence,  it  must  be  continuouSo   Similarly  the  second  term 


■  iimu 


leads  to  the  expression     l^  i'Lim}  C  -     j  j/Vil     M((H  t  Ol-~i)tj  16 


'  ^(h""]]' 


which  converges  to  a  continuous  functioHo 

The  last  term  in  (A)  is     ^^— ■  ^— "^  dt, 

Since   I  f  p/x)cJx|  "^  \  \^M   d/    ,  it  follows  that 


M 


.do 
e 

2nLu 


di 


=21^ 


.iiu 


r^  e Lm.^uy  ,.        /__L  1 

L  ^^'''        "^^       "^^     ^   K   mO[        u^j       by  corollary  b  to  theorem  10 

for  some  constant  ko 

Hence  this  term  approaches  zero  as  m-^oo   „  Therefore,  the  continuous 
functions  h^(u)  converge  uniformly  in  some  interval,  u.  ^  u  ^  u^,, 
to  a  continuous  function  h(u)o 

«(uj  -  H(u,)  =  Ur,    [     e ^ n,utj  Jt 


.lA  I  .  rm       -itu 


J   n^(uJJu      ^^'^^^  \(^)  converges  urdformly 
to  h(u)  the  limit  may  be  taken  inside  the  integral  and 

Since  H(u)  is  continuous  at  zero  and  we  have  shown  '  b solute  con- 
tinuity elsewhere J  it  is  absolutely  continuous  on  an  interval  in- 
cluding the  point  zeroo 

In  caseE|2,  j^oO  andir^2''+K'-l ,  integration  by  parts  of  (A)  and 
corollary  b  to  theorem  10  shows  that  the  first  r  derivatives  of 

h(u)  are  continuous  if  u  /  Oo  •  The  integral  expression  for  h(u) 

(r) 

in  terms  of  l|^  (it)  shows  that  ^>f-fl    implies  h^  '  (u)  is  continuous 

at  zerOo  This  implies  that  G(u)  =  a  ■l-J^^(vjd^   for  uJ-O;  g(v) 

is  continuous  for  v  ^  0,     If  EdjJ^^,  then  g^^^  (n)   is  continuous  for  u^O 

provided  r<3r+/<-|  and  is  continuous  for  u  s  0  provided  r<  y-|   „ 


.22- 


Corollary  to  theorem  11,  p   =  P(zn  s  r)-^0  uniformly  in  r,  r2  1 


nr 

as  n  -^  C30  . 


Since  G^^(u)™>G(u)  uniformly  for  u>0,  G(u)  is  uniformly  continuous 
for  u>0  and  G(u)  is  right  continuous  at  us  0^  the  corollary  is 
proveno 

It  is  interesting  to  note  that  Yaglom  [22]  treated  this  subject 
in  a  manner  similar  to  Harris o  Yaglom  defined 

^iLt]^  ]    6  a  n{\l]  and  proved  the  following  theoremo 

Theorem  l5o  li  ju>\  ,   the  variance  of  z,  is  not  zero,  and  f  (1)  4  oo 

then     rCil^iitJlj-ajfa]  ^(\-cl]  \^{^i]  ia 

-'''■■  In  this  paper  we  followed  the  general  line  of  Harris'  argument 
and  defined  k(s)  such  that  i^ljUS}^   KLL-Q-)  ^i^^s)  faj  -a   ^ 

This  led  us  to  the  conclusion  that  H(u)  was  the  distribution  function 
for  (l-a)w  which  is  the  limit  of  (l-Pj^Q)-{r 
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MEAN  NOT  GREATER  THAN  ONE 

We  have  examined  the  case  for  the  mean  greater  than  one  and 
now  we  turn  to  the  situation  which  corresponds  .to  the  problem  of 
subcritical  or  critical  atomic  pileSo  From  theorem  3,  the  prob- 
ability of  extinction  is  one  and  our  fissiordng  process  will 
eventually  die  out., 

When  the  mean  is  exactly  one,  the  limitirig  distribution  of  the 
conditional  random  variable,  y^.^,  (which  was  defined  in  the  last 
section  to  be  the  number  of  members  of  the  n^"  generation,  given 
that  this  number  is  not  zero)  is  the  same  for  any  discrete  branching 
processo 


DO 

F  (s)   =  "^  p"^_s^     be  the  generating  function  of  the  p'?! 

K  =   lim    l-iT^nlQj  J       if  the  limit  exists 

Theorem  16 „      If    /^  -  I      T^d)  ^(9  ,  f '"(')^oO  then 

(The  proof  may  be  found  in  Yaglom  [22] » ) 

Now  we  will  consider  the  limiting  distribution  of  Fj,^(s)   in 

the  case>a^lo 

-2U= 


Theorem  17 o  lfLL^\     and  f''<aO      then 


CO 


iir^  F^(s)=F(s) 


and   Fi:F(S)]=^F(s)+(l-yuj 

F(l)  =  /     f'(i)-K'-^ 


sU 


Proof: 


oo 


'^-   PTtjIT.^ ^^"^  ^y  definitiori  of  conditional 

rn  TK^n^O)  probability 

—  ^  P(Zr>  =  r)  r 

—  y   •  .-./  '^  -~  S  since  the  summation  is  over  r  0 


(1)    r  Is)  =  /  ,*.  k.(iH. 


The  following  theorem  may  be  found  in  Konig  [lU] o  If  Q{s)   is  analytic 
in  a  domain,  R^  such  that  @(s^)  -  s^  and  j  ^'(s^)  I  ~  z.C\   for  some 


I 

s-  in  R,  then,  in  this  domain   ^l£i--__~.£>_--^   converges  uniformly  to 

a 

a  continuous  function  D(s)    such  that       DLBiS)]  ■=  aVv>) 


f(s)   is  analytic  in  js|<:  !»      Since  f(l)   -  1  and  f'(l)<  1,    (/U<1), 
the  s^  of  Konig's  theorem  is  oneo     We  can  conclude  that       « r^(^''"_L 

converges  uniformly  to  a  continuous  function  D(s)    such  that 
^     (2)   D(f(s))   s/D(s),     D(l)   =  0  ,   d'(1)   -  lo 


^         From 


/      i 


(1).  F.,s)=i+(i^jri-)(v^ 


Lim      p  is)  ^F(s)  ="'^      (I    4 


jm 


('^(^--4(--s^)) 


(3)  F(s)  -  1  +  D(s)  K    by  Konig's  theorem,  since  convergence  of 


Lhh 


7^ 


n 


is  uniform,  the  limit  of     .1"^ 


exists c 


From  (3),  then,  D(s)  =  I  Fls)«|)|< 
From  (2),  ^(Pls))  =  ^  0(s) 


^^u/C[  F(s)-G 

F[f(5)]=  I  f  M^(li.'     by  (3) 


4  /^C 


KrF(s)-l] 


/< 


i  +  yC^  F  I  S  )  --y(/(. 


I 


And,   finally,       F  C  fls  1]  =-    >a  F  ( S  j  f  0  -/a) 
From  (3),     F(i)  -    I  ^    0(0  K"' 

0M1  =  £)      by(2) 


since 


f'(!)  =  4"FM 


o's 


IS^I 


4|l+    K"'0(SJ)J        :.yp'(S) 


ds 


fC"       by   (2) 


s- 
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One  of  the  more  interesting  questions  in  the   subject  of  branching 
processes  is  the  determination  of  the  number  of  generations  to 
extinction  whenyW^lo 

Let   ^i::  rnm  Irr^  2^,.  ^0] 


Theorem  1  8c  The  probabilities  b^^  satisfy  the  recursion  formula^ 


OS 

Proof  •  W     —  I  -  1"^ 


TheoreiT  I9.  The  moment  generating  function  of  N  is  given  by 


Proofs   e/s)  £  S'  I  p   -p  )  e 


:  /  -h  fe^»i)ej(5.) 

=27= 


So  ©(s)  is  determined  once  9-^is)   is  knoim.   If  we  let  z  "  e®, 
^-|_(s)  is  a  power  series  whose  coefficients  are  successive  iterates 
of  f^b)  =  1  =  f(l^b)  since  b^^^  s  f^(b^)  ^  f^.^Cb^),  where  f*(0)  =  0 
f^  (0)  =/^<  lo  (see  Fatou  [83o)  A  function  of  this  sort  is  meromor» 
pi-iic  with  poles  at  s  =  -  n  logyU  ^  n  s  1^2,.  „  o  So 

converges  everywhere  except  at  the  poleSo  (see  Lattes  [29]  ) 

Us)   s^  x^si  is  determined  by  X(,ixs)  -  f^Hx(s)) jX"(l)~  X'(l)=  1- 

y-Q   is  found  by  X(yo)  ^  b^  s  i  .,  p^  ©r  by  using  inverse  functions 

X='l(f*(s))  =//X"l(s)  and  determining  X^lCb^),, 

If  our  process  is  of  the  type  in  which  the  "parent"  does  not 

die,  then  the  total  number  of  individuals  is  of  interesto  If  m  i  1, 

the  probability  that  this  number  is  finite  is  equal  to  one  by 

theorem  iia^ 

Let  S'^  2^;   then  p(S<cP):s!   for  ^u^  I 

Let  q^,  =  P(5*r);  915)^2^^'"^°^  \sU^, 

Theorem  20.  Ifyu£  1^  then  g(s)  s  s  f(g(s))o 

r. 

Proof;  Let  S^"^  ]>  Itr^     It-^O    for  some  n  with  probability  one 
when  // ^.  1  by  theorem  ha.  Therefore,  Sj^  converges  to  some  random 
variable  S  with  probability  oneo  The  random  variable,  Sj^^  represents 
the  total  number  of  individuals  in  the  generations  from  the  zero^"  to 
the  n^h  generationo  Let  %^^   \   '^Vn   where  Yr-s^  "'?2r;   ?  since  I.,  =  5  o 
Let  6i^(Si  =^>Ci^,,S'''   .where  %c=P.   ^^^  1v>  =  PlVr^^'J- 
The  eveiit  y,^  -  n  can  be  divided  into  n  mutually  exclusive  events 
with  probabilities: 

P(z-L  s  i)P(n<=i  individuals  in  the  followirig  k-1  generations) 

for  is  1,2, o  o  c 
=28^ 


So  G(s)  is  determined  once  ©-i  (s)  is  knoirnic   If  we  let  z  =  e®, 
Q-^(s)   is  a  power  series  whose  coefficients  are  successive  iterates 
of  f-'(b)  =  1  ~  f(l==b)  since  b^^^  '   f*(b^)  =  f^-a^^o).  "^ere  f*(0)  =  0 
f*  (0)  -/A.<  i„  (see  Fatou  l83o)  A  function  of  this  sort  is  meromor- 
phic  with  poles  at  s  =  -  n  logyU  j  n  s  152^.  o  o  So 

converges  everywhere  except  at  the  poles^  (see  Lattes  [29]  ) 

X(s)  s^  x^si  is  determined  bj^  liu.s)  -   f^*(X(s)),x"(l)=  X"(l)^  1^ 

Yq   is  found  by  X(yo)  -  b^  s  l  ..,  p^  or  by  using  inverse  functions 

X-l(f*(s))  syur-^(s)   and  determning  X=l(b^)„ 

If  our  process  is  of  the  type  in  which  the  "parent"  does  not 

die,  then  the  total  number  of  individuals  is  of  interesto  If  >x  i  1, 

the  probability  that  this  number  is  finite  is  equal  to  one  by 

theorem  hsL^ 

Let    S-^^Zk,;   then    pfS^cQjsl      for     ^^  | 

Let     q^.  =  P(S^r).  9ls)"^^„/for      \sU\. 
■  '  tso 

Theorem  20 o    IfyU4  Ij   "then  g(s)   s  s  f(g(s))o 

_n 

Proof;  Let  S^  ~  ">  1h     2r<~0    for  some  n  with  probability  one 


L^i 


when  j}/ ^,   1   by  theorem  Ua„  Therefore,  Sj^  converges  to  some  random 
variable  S  with  probability  onCo  The  random  variable,  Sn>  represents 
the  total  number  of  individuals  in  the  generations  from  the  zero^^  to 


the  n^^   generationo  Let  S^^  *  "^Vn   where  VVs  ",>  2^^   ,  since  I 
Let  6«;S' =>q^,,S^''   ,  where  q^.^^Q^       and  5^.^=  P(yr  =  *^l« 
The  eveiit  y,^  -  n  can  be  divided  into  n  mutually  exclusive  events 
with  probabilities : 

P(z]_  ~  i)P(n=.i  individuals  in  the  followiiig  k»l  generations) 

for  is  1,2, o  o  c 
=28^ 


starting  with  one  individual  at  time  zero,  the  probability  of  pro- 
ducing n  individuals  in  the  next  k  generations  is  the  coefficient 
of  s^  in  Gj^(s)o 

Therefore,  q,    s      PCvt  e  i)p(n-i  in  the  next  k-1  generations) 

"  ■>  Pj_( coefficient  of  s^'=^  in  (Gj^^iCs))^) 

-  \  Pj_( coefficient  of  s^-  in  {s(}]^^]_(s))i) 

S.: 

'   coefficient  of  s"  in  (  ' ,.     Pj_(sGj^_3_(s))i) 
s  coefficient  of  s^  in  (  >  P^  (sGj^_i(s))i) 

s  coefficient  of  s^  in  f  (sGj^^2^(s)) 

Therefore^  Gj^(s)  =  '" ,"  <\ ^^   -  f  (5Gj^^3^(s))  o  Clearly^  gj^(s)  s  sGn(s) 

%  0 

since  P(yj^  =  n)  s  p(Sn  -  n-  1)  and  gj^(s)  "  sf (gj^^-^Cs)).  Taking  the 
limit  as  n  aporoache s  infinity^  we  see  that  g(s)  -  sf(g(s))o 

Usir^  theorem  20,  the  expected  total  number  of  members  of  a 
family  is:  ^'(0^4  cPf^'sll!   ^   h  [q/s^]  -f  S  ^  '  T  ^fsH^'iS)! 


'Ss!  ?S=J 


4  rhl^'l!) 


Then,    Q' U)  ^-  -r--—~- 


Similarly,  o'^fO  -■   L  '  '  ^ /^.VJ:^_l_   and  the  variance  of  S  isj 


^T   -:-   q'il)  T  Q^i:}  ■•«  [9''''] 
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ESTIMATION 
In  the  introduction  we  mentioned  that  the  parameters,  p^^   must 
be  empiraca  ly   determined  for  the  given  amount  and  arrangement  of 
atomic  material.  The  method  of  maximum  likelihood  can  be  applied; 
however,  it  seems  necessary  to  assume  that  certain  random  variables 
are  observable. 

Let  Zj^jp  be  the  number  of  individuals  in  the  m^^  generation 
who  have  exactly  r  descendents  in  the  (m4l)^^  generation.  Then 

2m -^  ^mr  ^   ^   "^^  Let  S^  ^^  1^ 

Theorem  21,  The  maximum  likelihood  estimates  of  p^*  and  jU    based  on 
observed  values  of  z^^p  for  m£n  (n-f  1  observed  generations)  are 

Proof:   In  order  to  construct  the  likelihood  function,  first  we 
determine  the  conditional  distribution  of  Zj^,  r  .  0,1,  ,  o  ,  given 
Zjjj  o  From  our  initial  assumptions,  each  member  of  the  m^*^  generation 
has  the  same  probability  distribution  associated  with  his  progeny. 
Therefore  the  probability  of  k^  members  having  no  descendents  is 

P  ^  I  k*^/  "^^^^^     jT'  1  ^^  ^^®  binomial  coefficient  and  is  equal  to 

-jzi rrVnT"!  -  Of  "t^he  remaining  Zr«  -  k^  members,  say  z^  -   k-.  have 

one  descendent.  This  probability  must  be  |  "^   J  P    ,  Continuing 

in  this  maruier  we  get  the  multinomial  functions: 

P(2mo3  Zral'  °  ■•  »>  2^>  '  °  °|  ^m-l,0>  2^„i  ^  ,  o  o  .,  z^_-^  ^,    <, .  „) 

7     J  7T  p^'M-K 
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^- 


Therefore  P(z^o,  z^i,  .  .  .  |  Zn.1,0.  ^n-l,!.  =  °  ')^^h,-l,0>    ^^.i^i 


2 


lr-0 

Due  to  the  Markov  properties  of  the  branching  process  this  is  the 

joint  distribution  of  the  Zj^  ,  m  =  0,1,.  ,  .  ,  r=  0,1,  .  ,  . 

The  only  factor  of  the  likelihood  function  which  depends  on  the  p^,  is 

oo    >  7 
If    ?r 


oo    n_  \ 

the  logarithm  of  this  expression  is  ^  |  ^  Z^^    ^'^fe^fVv 

oo 
To  maximize  this  expression  subject  to  the  restraint,  ^'    P  =  j 

the  method  of  Lagrangian  multipliers  may  be  used. 

Setting  AL.  —  0  ,  we  get  ^  ^    —  Jp  •^^"'^  r  -  0,1^  o  o  . 
But  ^"  P  -  0  ',   summing  the  above  expression  over  r  gives: 

The  double  sum  of  the  Zj^  is  merely  Sj^  so  <^  =  S^^  and 

h 
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CO  »« 

To  find  Ju    recall  that  /^  =  2^  rp^  and  that  ^  l^Zj^^  =:  Zk^,+  ( 
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Using  these  facts,  it  follows  that: 
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ION  FLUCTUATION 
The  number  of  ion  pairs  produced  by  a  fast  primary  particle 
passing  through  an  absorbing  medium  is  a  branching  stochastic  process 
(Moyal  [18]  ).  In  this  process  the  primary  particle,  which  may  be 
a  gamma  ray,  has  a  certain  charge,  mass,  and  velocity  or  a  given 
frequency J  the  thickness  and  atomic  properties  of  the  absorbing 
medium  are  known.  The  primary  interacts  with  the  atoms  in  the 
absorber,  freeing  electrons  by  the  Compton  effecto  These  electrons 
are  the  first  generation  and  each  may  ionize  another  atom  or  atoms 
depending  on  their  energy;  secondary  electrons  may  ionize  more 
atoms,  etc.  The  process  ends  when  all  the  freed  electrons  are  too 
slow  to  produce  further  ionizationo  Tiiis  results  in  some  number  of 
free  electrons  and,  thereby,  an  equal  number  of  ion  pairs. 

In  determining  the  distribution  of  this  number  we  shall  make 
the  following  assumptions: 

1„  Successive  ionizing  collisions  are  statistically  independent. 

2„  The  total  energy  loss  of  the  primary  is  very  much  smaller 
than  its  initial  energy  and,  hence,  the  loss  of  primary  energy  may 
be  neglected.  Therefore,  we  are  not  treatitig  slow  primaries  or 
very  thick  absorbers.  We  also  assume  that  energy  losses  due  to 
radiation,  nuclear  interactions,  etCc  are  negligibleo 

3o  The  absorbing  medium  is  homogeneous. 

Uo  The  recording  of  data  is  delayed  sufficiently  to  allow 
tenaination  of  the  process.  Photoelectric  effects  in  the  gas  or 
chamber  walls  is  neglected. 

Unlike  the  branching  processes  treated  in  the  previous  sections 


the  probability  distribution  of  the  number  of  progeny  per  individual 
in  a  given  generation  is  dependent  upon  the  generatioHo  Thus  one 
may  expect  different  methods  to  be  used^  Moyal  has  applied  the 
Laplace  transform  method  to  solve  the  Kolmogorov  equations  of  a 
Markov  process  and  this  is  outlined  beloWj 
Let  t  be  the  thickness  of  the  absorber 

N  =  the  number  of  absorber  atoms  per  unit  volume 

6"(E)  =  total  cross  section 

q  =  primary  ionizing  collision  rate 

qj.^  =  P(n  ion  pairs  are  produced  in  a  given  collision) 

p  (t)  =  P(  n  ion  pairs  produced  in  a  thickness  t) 

__  >■■' 

M(s,t)  =  "^  pjj(t)e'"^^^  ,  the  generating  function  of  Pi^(t) 

Then  (n  6'(£)  (dl^  dt)  is  the  probability  of  an  energy  loss  between 
E  4nd  £  +  dE  in  an  absorber  of  thickness  dto  The  primary  ionization 
rate  is     Cj  =  (\j  \  (T  (^1  dE 

V 

The  ^^  9(0  ~  ^   ^"^  q(dt)  =  F(  an  ionizing  collision  in  thickness  dt) 
Since  at  least  one  ion  pair  is  produced  in  a  given  collision,  q^  =  Oc 
Based  upon  the  above  assumptions  we  have 

W  PJt,+tj)  =^  P(t,)p  Itj.) 

Kronecker  delta  and  0(St)   goes  to  zero  faster  than  gt  as  Jt-^0. 
By  assumption  1  and  equation  (A), 


"^[^0   i-(M(5.St)-ljM(s,t) 


Comparing  this  equation  with     .PX.^.  —  cy   ,   it  is  clear  that 


oo 


M 


QRls) 


where      Q^C^t  ^^^)  ~  ^  {e''^^''\)% 

The  inversion     of  the  Laplace  transfomij,  M(s^t),   is  given  by: 
Pfl\Q)  ~  Ywl     ]  ^  e        do  (see  Widder  [2l] ) 

Evaluating  this  integral  by  the   saddlepoint  method  (see  Korn  [l^]), 
using  the  first  term  only,   gives: 

Pn(Ql  =4-[;jTTQR"(s„)]"''''e>fp[Q(Rls«)-S„R'(S„)] 

where  S^  is  related  to  n  by 
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Since   R(S„^  =^  (e    -l)^K 


and  c  is  a  normalizing  constant. 
Since  ^  p^(Q)  =  I 
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